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ABSTRACT 

Galaxy evolution is driven to a large extent by interactions and mergers with other 
galaxies and the gas in galaxies is extremely sensitive to the interactions. One method 
to measure such interactions uses the quantified morphology of galaxy images. Well- 
established parameters are Concentration, Asymmetry, Smoothness, Gini, and M20 of 
a galaxy image. Thus far, the application of this technique has mostly been restricted 
to restframe ultra-violet and optical images. However, with the new radio observatories 
being commissioned (MeerKAT, ASKAP, EVLA, WSRT/APERTIF, and ultimately 
SKA) , a new window on the neutral atomic hydrogen gas (H i) morphology of a large 
numbers of galaxies will open up. The quantified morphology of gas disks of spirals 
can be an alternative indicator of the level and frequency of interaction. The Hi in 
galaxies is typically spatially more extended and more sensitive to low-mass or weak 
interactions. 

In this paper, we explore six morphological parameters calculated over the extent 
of the stellar (optical) disk and the extent of the gas disk for a range of wavelengths 
spanning UV, Optical, Near- and Far-Infrared and 21 cm (Hi) of 28 galaxies from 
The Hi Nearby Galaxy Survey (THINGS). Though the THINGS sample is small and 
contains only a single ongoing interaction, it spans both non-interacting and post- 
interacting galaxies with a wealth of multi-wavelength data. We find that the choice 
of area for the computation of the morphological parameters is less of an issue than 
the wavelength at which they are measured. The signal of interaction is as good in 
the Hi as in any of the other wavelengths in which morphology has been used to 
trace the interaction rate to date, mostly star-formation dominated ones (near- and 
far- ultraviolet). The Asymmetry and M20 parameters are the ones which show the 
most promise as tracers of interaction in 21 cm line observations. 

Key words: galaxies: fundamental parameters galaxies; spiral galaxies: structure 
galaxies: interactions galaxies: kinematics and dynamics 



1 INTRODUCTION 

Evolution of galaxies in the cold dark matter with a cos- 
mological constant Universe (ACDM) appears to be driven 
by the merger and interaction of dark ma tter haloes (e.g., 
the Millennium Simulation by ISpringel et a l. 2005). There- 
fore, a substantial observational effort has been made to 
quantify the rate of mergers and interactions over time. Sev- 
eral methods have been developed to estimate the interac- 
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tion rate: identification of physically close pairs of galaxies 
in re dshift surveys (e.g, iPatton et al.l I2OO0I : Ide Ravel et al.l 
l2009l ). measures of H-alpha equivalent width, far-IR flux 
from (ultra) luminous infrared galaxies ((U)LIRGs, see 
e.g., iMurphy et al.l [200lh. frorn the s tar- burst, OH mega- 
masers (e.g., Klockner fc BaanI l2005l : iDarling fc Giovanellil 
l2006h and ident ification of galaxies with d i sturb e d mor- 
phologies (e.g. IConseUcel l2003l : iLotz et al.1 |2004 l2008al : 
IConsehce et al.l |2009| )." This observational effort has been 
matched by theoretical ones to accurately map the well- 
understood merger trees of Dark Matter haloes onto galaxy- 
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galaxy merger rates (see the review in lHopkins et al.]|2010l . 
and references therein). Thusfar, theoretical models suffer 
from large systematics. However, ongoing efforts in both cos- 
mological hydrodynamical simulations and semi-analytical 
models can be expected to match the observational accu- 
racy soon. 

Merger rates from disturbed morphologies of galaxies 
have been explored extensively with quantified classifica- 
tion of galaxies. Based on a series of scale-invariant param- 
eters, quantified galaxy morphology has been applied pre- 
dominantly on restframe ultra-violet images of galaxies (e.g., 
iTavlor-Mager et al.|[2007h . The advantages for optical or ul- 
traviolet are that interacting galaxies are star-forming and 
hence bright in ultra-violet and blue side of the optical. Their 
morphology shows enhanced surface brightness and clear 
signs of disturbance. The high surface brightness ensures 
more complete samples for a given observation. Observation- 
ally this approach is also attractive as it has the advantage 
of similar spatial resolution at high- and low-redshift, using 
the Hubble Space Telescope and Galaxy Evolution Explorer 
(GALEX) or Sloan Digital Sky Survey (SDSS) respectively. 
Disadvantages of this method are the time-lag for the in- 
teraction to trigger star-formation, and modification of the 
morphology due to dust obscuration. 

However, new windows for quantified morphology will 
be opening up; the far-infrared emission from star-formation 
and molecular gas are now resolved with Herschel and 
in the near future with the Atacama Large Millime- 
ter/submillimeter Array (ALMA). The atomic gas with 
its many signatures of interaction (tails, bridges, beards, 
warps, clouds etc.) will also be much better resolved with 
the commissioning of th e Square Kilometer Array (SKA; 
ICarilh fc Rawlingsi 120041 ') . and its precursors , the South 
Afric a n Kar o o Ar r ay Telescope ( MeerK AT; I Booth et al.l 
I2OO9I : iJonad I2OO7I : Ide Blok et al] |2009|). the Australian 



SKA Pathfinder (ASKAP: I Johnston et all bOOSl ). and_the 
pathfi nders, the Extended Very Large Array fEVLA: lNapieil 
I2OO6I ) , and the APERture T ile In Focus instrument (APER- 
TIF; IVerheiien et al.l |2008l ) on the Westerbork Synthe- 
sis Radio Telescope (WSRT). Spatial resolutions of these 
observations will start to rival those in the ultra-violet. 
Hi morphology also appears to be very sensitive to the 
smaller interactions with tidal features often reported to be 
much more visible than in any other waveleng th (see the 
"Hi Rogues Gallery'^ compilation in Hibbard et al. 200 L). 
Hence the morphology of galaxies in other wavelengths, es- 
pecially Hi might be an equally good or surpassing indi- 
cator of tidal interaction than that of the ultra-violet and 
other star-formation dominated emission. Notably, the en- 
visaged all-sky surveys with ASKAP (WALLAB'vH), and 
the WSRT/APERTIF will then provide an accurate census 
of mergers in the local Universe. 

The aims for this paper are to explore (I) which 
wavelength shows the strongest signal of interaction, (2) 
which morphological parameters are the optimal discrim- 
inators for interaction, and (3) over which area morpho- 
logical para meters need to be computed. In a companion 
letter (Holwerda et al.ll2010l 'l. we briefly highlight how well 
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Hi morphology shows the signal of interaction compared 
to the UV, optical and FIR. In further papers, we de- 
flne the Hi parameter space to identify interacting galaxies 
(|Holwerda et aLllioilbl ). using a sub-sample of WHISlH, we 
derive a ti me-scale for interactions to reside in this param- 
eter space (|Holwerda et al.ll2011cl ). and we infer the first 
interaction fraction an d rate based on the WHISP survey 
([Holwerda et a l.' 2011a). The Hi morphological phenomena 
in the Virgo cluster environment (e.g., ram- pressure strip- 
ping) are explored in iHolwerda eraLr(|2011dh . 

In ^we discuss our sample and data used, in section 
[3] we discuss the morphological parameters, as well as ef- 
fects of uncertainty and possible biases. Our results are in 
^ together with notes on each individual galaxy. Our dis- 
cussion on the suitability of the H I parameters is in and 
our conclusions are summarized in ijT] 



2 SAMPLE & DATA 

We use the publi c datasets from the H I Nearby Galaxy 
Survey (THINGS: [Walter et al . 20 08f|, the Spitzer Infrar ed 
Nearby Galaxies Survey (SINGS, iKennicutt et all l2003l fl. 
the GALEX Nearby Galaxy Atlas (NGAf| and the Sloan 
Digital Sky Survey Data Release 7 (SDSS, DR7fl Table 
A2 {in the online version) lists the availability of Spitzer, 
GALEX a nd SDSS data for ou r sample, as well as the basic 
data from IWalter et al.l (|2008l ) . Optical data is either from 
the SINGS ancillary data, Nasa Extragalactic Databas^fl 
public data o r SDSS, preferring SPSS. W e ai m to cover the 
sample from T rachternach et all (|2008') and |de Blok et al.l 
(2008) because for these galaxies there is detailed dynamical 
information (rotation curves, dynamical centre, and param- 
eterizations of non-circular motions). We require Hi data 
from the THINGS sample, Spitzer data from the SINGS 
sample of a program with equivalent quality data, and at 
least some op tical data. These re quirements narrow the 34 
galaxies from IWaher et all (|2008l ) down to 28. 

The H I column density maps are the naturally-weighted 
(NA) and robust-weighted (RO) total intensity maps, ex- 
pressed as Hi column density using the formalism and 
beam sizes reported in I Walter et al. I (|2008l ). We chose the 
naturally-weighted maps to define our contours as these are 
more sensitive to larger scale structures (of the order of the 
disk) than the robust weighted maps. The natural-weighted 
maps have lower spatial resolution (typically of the order of 
10" resolution) than the robust- weighted ones (typically 6"). 
We compute morphological parameters for both Hi maps. 
To define the area over which parameters are computed, we 
picked two Hi column density contours: 0.3 and 30 x lO^'' 
cm^^. These correspond to approximately the spatial ex- 
tent of the H I and stellar disk respectively, but may exclude 
areas corresponding to "Hi holes". 

The majority of the Spitzer IRAC and MIPS data are 
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from the SINGS project with t hree additions fr om the Local 
Volume Legacy Survey (LVL. iLee et"aLl |2008| ) : NGC 5236 
(M83), NGC 5457 (MlOl) and IC 2574. Spitzer data includes 
all IRAC (3.6, 4.5, 5.6 and 8.0 fim) and MIPS (24, 70 and 160 
yitm) channels. These tend to comfortably cover the stellar 
disk and most of the Hi outer contour (3 x 10^® cm"'^). 

In the case of the SDSS optical data, we obtained the 
original tiles around our galaxies and combined them into 
larger mosaics using swarp (http://astromatic.iap.fr/). 
The images were sky-subtracted before combination to ac- 
count for the different sky values in each scan-strip. The 
SDSS mosaics have distinct advantages over the NED and 
SINGS ancillary data: uniform depth, a well-defined set of 
filters and a field-of-view that covers the whole H I map. 

Most of the GALEX data is from the Nearby Galaxy 
Atlas (NGA) supplemented with the all-sky survey in two 
cases: NGC 3184 and NGC 6946. NGA data generally means 
the galaxy is in the focus of GALEX with 4" resolution in 
FUV (1528 A) and NUV (2271 A). However, galaxies such 
as NGC 3031 and M81 Dwarf A are part of the M81/M82 
group portrait and M83 is in the corner of the GALEX FOV 
because of a bright foreground star and these galaxies are 
therefore slightly out-of-focus with GALEX with resolutions 
closer to 6" 

Table Al {online edition) lists the resolutions and wave- 
length of all our data. GALEX resolution ranges from 4" -6" 
depending on position in the field, MIPS at 24 fim is 6" and 
the RO maps are typically also 6" resolution, depending on 
position in the sky and axis. The MIPS 70 and 160 ^im and 
the NA maps are of poorer resolution. 



2.1 Data Preparation 

To start, we shift the natural weighted Hi map, such that 
the centre of the galaxy is in the centre of the image. Because 
the THINGS observations are pointed, the galaxy is already 
close to the central part of the image and the shifts are 
sr nall (a few pixels). W e use the central positions reported 
in I Walter et al.l (l2008f). who lis t the d ynamical central posi- 
tions from Trachternach et al] l|2008l ). supplemented on oc- 
casion with Spitzer 3.6 fim central positions. Subsequently, 
we align all the different data using wcsmap and geotran in 
IRAF to this centered natural weighted (NA) Hi map. Our 
next step is to convolve the optical and IRAC data to 6", 
approximately the resolution of the majority of the rest of 
the data, before determining the morphological parameters. 
After alignment and smoothing, a mask of the 

foreground stars is cr eated using a Source Extractor 

l|Bertin fc Arnoutdl 19961 : lHolwerdall2005l ) catalog and result- 
ing segmentation map of the SDSS-i image, selecting small 
objects in the field. If the SDSS-i image is not available, we 
use the IRAC channel 2 (4.5 /^m). Channels 1 and 2 of the 
IRAC instrument trace the old red stellar component and 
we chose channel 2 as it does not contain the 3.1 /im PAH 
feature so any hot ISM region that belongs to the galaxy is 
spuriously rejected. The disadvantage of the IRAC data is 
that it may not cover the entire gas disk but to mask fore- 
ground stars, it is preferable to use optical or near-infrared 
rather than GALEX data. 



3 MORPHOLOGICAL PARAMETERS 

The morphological parameters we use here have been es- 
tablished with repeated applications to deep Hubble Space 
Telescope images and reference local galaxy samples. There 
is the Concentration- Asymmetry-Smoothness classification 
scheme and the Gini and M20 parameters. Ell ipticity is 
somet i mes added. We use the d efinitions shared bv | Conselicel 
(l2003h . lLotz et al.l (|2004l ). and lScarlata et all (|2007l ) for our 
computation of the six morphological parameters: Concen- 
tration, Asymmetry, Smoothness, Gini, the moment of light 
(M20) and Ellipticity. 



3.1 CAS 



[Abraham etHI (Il994 Il996al lbl) introduced definitions of 
asymmetry, concentration and contrast to classify galax- 
ies in the Hubble D eep F ield N orth. Followin g the w ork 



bv iBershadv et al 



200(t and iTruiillo et all (|200ld lbh. 



IConselice et al. ( 200C ) refined these parameters, culminat- 



mg m iConsehcd (|2003l ) which added a local volume refer- 
ence in R-band (the sample from iFrei et akl 1 1996!). The 
thus established parameter space has been used exten- 
siv ely on all HST wid e an d deep surveys; e.g, i n GO ODS 
bv lBundv et all ||2005| ) and iRavindranath et a.1.1 (I2OO6I). the 
HUDF by dYan et all l2005h. C OSMOS bv IScarlata et al] 



|2003), GEMS bv [jogee et all llioogil . and the extended 



Groth strip bv ITruiillo et al.1 1|2007^ !^ 



3.1.1 Concentration 

Concentration is defined bv ITruiillo et alj (|2001al ) as: 



C = 5 



Vr2o / 



(1) 



with Tf is the radius containing a percentage / of the light of 
the galaxy, in this case 80 and 20 percent respectively. Other 
definitions using rgo and r^o are also in vogue, notably in 
the S DSS catalogue (see the discussion in lCraham fc Driveij 
l2005h. The radii are often taken from Source Extractor 
llBertin fc Arnoutd Il996l : iHolwerdal l2005h output but we 
computed these here without the use of the Source Extrac- 
tor program. The two apertures for these radii are circular 
and hence this parameter is somewha t sensitive to the in- 
clination of the s piral disk (See also IScarlata et al ] I2OO7I : 
iBen do et al ] |2007l . and i ]3.5|l . Concentration depends on the 
adopted central positionfor the measurement apertures. 



3.1.2 Asymmetry 

Following earlier work bv lAbraham et al.l (Il996bh . the now 
most commonly used definition of Asymmetry is from 



most commonly 
IConselicd (|2003l ): 



A 



S».,-|J(i,j)-Ji8o(i,j)| 
S.,,|/(i,j)l 



(2) 



where I{i,j) is the value of the pixel at the position i,j in 
the image, and /igo(i, j) is the value of the pixel in the same 
position in an image, which is rotated 180° around the cen- 
tre of the galaxy. To compute Asymmetry, we need a known 
position of the centre of the galaxy as well as well-defined 
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area. lAbraham et"ai] (|l996bl ): [Coriselicel l|2003l ) apply a fur- 
ther correction to remove a contribution to the Asymmetry 
value by the sky background. See section [3.4.1l why we chose 
not to in this paper. 

3.1.3 Smoothness 

Following iTakamival ()l999l l. IConselicd (|2003l ) introduced 
Smoothness, which has gone through several definitions. 
Here we use: 

_ ^.j\I{i,j)-Is{iJ)\ 

where Is{i,j) is the same pixel in a smoothed image. 
Smoothness is a parameter ized version of the unsharp mask- 
ing technique iMahnI | 19781 ) used on photographic plates to 
identify faint structures. The various definitions employ dif- 
ferent smoothing kernels and sizes, the most recent one using 
a flexible kernel-size of 0.2 Petrosian radius and the boxcar 
shape. To simplify matters, we cho se to use a fixed 6 " Gaus - 
sian smoothin g for our definition. lAbraham et all (|l996bl ): 
IConselicd ((2003,) apply a further correction for a background 
contribution to this parameter, which we do not (See section 
13. 4. li as to why). 

3.2 Gini and M20 

The Gini and M20 parameters were established by 
iLotz et al.l l|2004l ) as an alternative to the CAS space. 



the position of the galaxy's centre. The total second order 
moment of an image is defined as: 

Mtot = SAf, = E/4(a;, - Xcf + {y^ - ycf]. (5) 

When we now rank the pixels by value, we can define the 
relative second order moment of the brightest 20% of the 
flux: 

M20 = log , for which S*/, < 0.2 Itot is true. (6) 

y Mtot J 

Pixel k is the pixel marking the 20% point in the list of 
pixels ranked by flux value. 

Some authors var y the central pos i tion (xc, Vc) to min - 
imize this parameter IjLotz et al. I I2OO4I : iBendo et al.ll2007l '). 
Because we have dynamical centres, we flx Xc,yc and treat 
deviations from this value as a source of uncertainty. 

3.2.3 Ellipticity 

IScarlata et al.l (|2007l ) added the ellipticity of a galaxy's im- 
age to the mix of parameters in order to classify galaxies 
according to type in the COSMOS field. Ellipticity is de- 
fined as: 

E = l-h/a (7) 

with a and h, the major and minor axes of the galaxy, com- 
puted from the spatial second order moments of the light 
along the x and y axes of the image in the same manner a s 
Source Extractor l|Bertin fc Arnouts|[l99^ : iHolwerdalbOOSl ). 



3.2.1 Gtm 

The Gini parameter is an established qualifier in economics 
for the inequity in income for a population. For a Gini value 
of one, every person (or pixel) owns an equal fraction of the 
wealth (or fiux). A Gini value of zero, all wealth (or fiux) 
is con centrated in a single person (or pixel)^ Abraham et al.l 
l|2003t ) introduced the Gini parameter, and lLotz et al.l ( 20041 ) 
used this scale-invariant parameter to characterize the ho- 
mogeneity of a galaxy image. This parameter shares some 
of the characteristics of Concentration and Smoothness from 
the CAS space but does not depend on the size and shape 
of a convolution kernel or the choice of the galaxy's centre. 
In ILoIz ct al. (200j|), their equation 6, following the work 
bv iGlasseii (|1962l ). the Gini parameter can be redefined for 
speed if one orders the pixels according to value first: 



G 



In(n 



-Ei(2i-n - 1)1/^1 



(4) 



where li is the value of pixel i in the ordered list, n is the 
number of pixels in the galaxy image, and I is the mean 
pixel value in the image. We implemented this definition of 
the Gini parameter as the computationally least costly one. 



3.2.2 M20 

iLotz et ai] l|2004l ) also introduce the relative second order 
moment of the image to classify galaxies. The second order 
moment of a pixel is: Mi = Ii[{xi — XcY + {Vi — i/c)^], where 
li is the value of pixel i in the image, and Xi and yi are 
the X and y coordinates of that pixel and Xc and yc are 



3.3 Uncertainty Estimates 

The sources of uncertainty in the above parameters are: (1) 
shot noise in the pixel values, (2) uncertainties in the po- 
sition of the centre of the galaxy and (3) variations in the 
area over which the parameters are computed. The first two 
uncertainties can be estimated using a Monte-Carlo set of 
iterations, the last one using a jackknifing technique. The 
relative contribution of these sources of uncertainty depends 
on the instrument characteristics and hence wavelength, res- 
olution and distance of the object. For instance, photon shot 
noise is more of an effect in UV and optical data, compared 
to the Hi. Here however, we treat each image the same and 
compute these uncertainties for the relevant parameters. 

The shot noise effect on a parameter can be estimated 
by reassigning random pixel values around the mean value 
to each pixel in the image and recomputing the parameter 
several times. With a few iterations, the rms of the spread 
in parameter values is an estimate of uncertainty in the pa- 
rameters. 

Similarly, the uncertainty due to the measurement error 
of the centre of the galaxy {xc,yc) is can be estimated by 
varying the input centre, recomputing the parameter and 
calculating the spread in values over a certain number of 
iterations (ten in our case). We use random deviations from 
Xc and yc within a normal distribution with a width of 6" 
to mimic the uncertainty in the position of the centre. 

The latter uncertainty estimate is important for Con- 
centration, Asymmetry, M20, and Ellipticity as these depend 
on the assumed centre of the galaxy. To minimize our dif- 
ferences with dynamical parameters later, we adopted the 
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dynamical centre from iTrachternach et alj l|2008h whenever 
possible. In previous optical high- ft studies, the uncertainty 
in the position of the centre was considered less of an issue 
-since it is generally much less than a resolution element- 
but we find that for nearby galaxies, this is a substantial 
part of the error-budget in the parameters. 

The Gini parameter does not depend on the central po- 
sition and it is only weakly sensitive to shot noise in the 
pixel values. Therefore, we computed its uncertainty from a 
shot noise Monte-Carlo set and the rms in Gini values from 
a series of subsets of the pixels in the image (a jackknife set). 
The jacknife approach t o the Gini parameter uncertainty is 
advocated by lYitzhakil |l993). Gini shot noise and jacknife 
uncertainty estimates are of similar magnitude. 

The reported uncertainties in Tables (A3-A30, online 
edition) are the combined uncertainty from shot noise and 
central position for all parameters, with the exception for 
the Gini parameter. In the case of Gini, the uncertainty 
is the combination of shot noise estimate and the jacknife 
estimate. These are formal errors and in our opinion, the 
actual errors in these parameters are larger, predominantly 
due to viewing angle, resolution effects, image artifacts, and 
remaining uncertainty about which pixels to include in the 
computation. 

3.4 Systematics 

There are several effects that may influence the structural 
parameters: (1) the sky background, notably any structure 
in this background, (2) signal-to-noise, (3) the inclination of 
the disk and (4) the resolution and sampling of the instru- 
ment or conversely the distance of the galaxy. These system- 
atics vary in prominence for each different wavelength but 
we focus on the effects on Hi morphology here. 

3.4-1 Background Contribution 

The sky and instrumental background in these images in- 
fluences the morphological parameters in three different in- 
stances. First, the choice of area over which the parameters 
are computed is influences if a s/n criterion is handled for the 
inclusion of pixels. Secondly, the mean sky background over 
the area contributes to the total parameter value. This is a 
de-facto weighting of the shape contribution to the parame- 
ter. And thirdly, in the case of Asymmetry and Smoothness, 
the noise in the background adds value since these parame- 
ters use the absolute difference in pixel values. 

The first consideration does not come into play in our 
experiment as we fixed the choice of area based on the H I 
contour. The contribution by the mean background is min- 
imized in our case as we selected uniform data, which is 
already sky-subtracted, and we carefully sky-subtracted the 
SDSS tiles before addition. Especially the SINGS data was 
meticulously corrected for backgro und contributions ( see the 
fifth SINGS data release notes, ISINGS teaiiil |2006| ). The 
mean background contribution is dominated by the shape 
over which the parameters are computed. To compare, we 
include the parameters computed for the shape with all pix- 
els set to unity in Appendix A (online version of the paper, 
Figure A1-A27). This is effectively a gross overestimate of 
a mean background contribution. This leaves the noise term 
in Asymmetry and Smoothness as a separate problem. 



iLotz et al] l|2004 ) correct their Asymmetry and 
Smoothness values for a background contribution. Because 
these parameters use a absolute difference between pixels, 
any background noise -or other background contribution- 
will add to the signal in both these parameters. To correct 
for this, most authors estimate the background contribution 
in an area of similar size as the object near or around the 
object. In the case of Hi maps, this would be complicated 
by the fact that the H I map covers most of the field of view 
in the other wavelengths. 

The reason we did not choose to correct our parame- 
ters for a background contribution this way was because the 
annulus where the background contribution would be com- 
puted for the inner contour would fall within the outer H I 
contour. Hence, the background may contain signal we are 
trying to measure. Using an annulus further out of the tar- 
get galaxy is complicated by a lack of coverage in the optical 
and Spitzer images. 

To test our assertion that the background contribution, 
which includes any fiux by a background source that was 
not masked, is minimized, we compare the parameters in 
both contours in IRACl (3.6 fim) and MIPS 24 ^im im- 
ages, two filters for which we have measurements for all our 
galaxies. The 3.6 images were smoothed to 6" resolution, 
correlating the noise which would strengthen any weak back- 
ground patterns and the MIPS 24 was not smoothed but a 
likely candidate for background structure. Figure [T] shows 
the comparison between the two contours. While there is 
much scatter, there is no systematic offset in either Asym- 
metry or Smoothness. Similar plots for the other wavelength 
show a similar trend. This represents a lower limit to a back- 
ground structure contribution from for example instrument 
characteristics etc. However, the lack of an offset is encour- 
aging. The shot noise uncertainty ( ^3.3|l is separate effect 
from any background contribution but we note that, espe- 
cially in Asymmetry, any background deviations seems to 
fall within the computed uncertainty. 

Alternatively, one could construct background images, 
mimicking the background mean value and noise. However, 
if these do not include a mean background contribution, 
there is very little change in the values for all morphological 
parameters (at least for the 3.6 and 24 /.tmimages where we 
tried this) . The large signal-to-noise over many pixels in the 
selected shape dominates the parameter value. 

Therefore, we do not correct our values for a background 
contribution here. However, these parameters have been de- 
termined for a large number of pixels and high signal-to- 
noise images of nearby galaxies. As soon as the background 
can be a substantial contribution to the intensities in an 
image, such is the case in HST images of distant galaxies, 
the background contribution needs to be revisited. In the 
case of H I observations in particular, this issue will be much 
less prominent as the column density map is constructed 
from high signal-to-noise line detections, lowering the rela- 
tive contribution of the background. 

3.4.2 Signal-to-Noise 

iLotz et al] (|2004l ) also conduct an experiment how much 
these parameters are affected by a change in signal-to-noise. 
They find that G, M20, and C are reliab le to w i thin 1 % for 
galaxy images with S /N > 2 per pixel. iLiskeij l|2008l ) find a 
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Figure 1. The Morphological parameters for the smoothed IRACl (3.6 lira) and native resolution MIPS (24 /xm) images for the gas 
disk contour and stellar disk contour. If there is a significant contribution by a background noise to any of the parameters, this should 
manifest itself as an offset between the two contours. 



similar signal-to-noise limit for the Gini parameter. Asym- 
metry and Smoothness decrease with S/N but stay within 
Umits above S/N > 5. The stellar disk of all our galaxies are 
detected well beyond this noise level. The Hi perspective 
will initially be used in the nearby universe where the signal- 
to-noise level will comfortably exceed these limits because 
line detections exceed s/n of 3 or 5 in individual channel 
maps to be included in the column density map. Signal-to- 
noise effects will need to be considered when the classifica- 
tion scheme is applied to galaxies at higher redshifts with 
for instance SKA. The redistribution of pixel values around 
the mean (the Monte-Carlo error estimate) in our parame- 
ter is in part to quantify how much effect noise has on our 
parameter determinations. 

3.4-3 Distance Effects on Morphology 

The accuracy of the morphological parameters depends on 
sampling and resolution of the analysed images. In addi- 
tion, cosmological surface brightness dimming plays a role 
at greater distances. Ifere, we focus on the effects of dis- 
tance on the morphological parameters measured in our H I 
column density maps. In this section, we focus on the ef- 
fects on increasing distance on the morphological measure- 
ments in our Hi column density maps, to ensure that these 
THINGS measures do not need to be corrected for relative 
distance. 

Distance effects in the optical and inf rared images have 
been di scussed by b oth Lotz ct alj (|2004h and iBendo et al.] 
((2003). [ Lotz et al] |2004) find that decreased sampling 
has the strongest effects on Concentration and M20; 15% 
changes when pixel scales become greater than 500 pc. Gini, 
Asymmetry, and Smoothness, on the other hand, are rela- 
tively stable with decreasing spatial scales down to 1000 pc. 



In the near and far-infrared, iBendo et al.l (|2007l ) find that 
Concentration, Gini and near-infrared Asymmetry (yla.snm.) 
are invariant with image smoothing and the M20 parameter 
is only moderately affected (<20%). The one exception is 
the Asymmetry in the far-infrared (24 /xm). They report a 
dramatic change with distance; A24,^m. oc with d in 

Mpc. The spatial resolution of these 24 pm images is lower 
than the images at other wavelengths. However, they do 
sa mple smalle r phys ical scales than the limit 1 kpc found 
bv lLotz et al.l |2004) for the applicability of these parame- 
ters. The dramatic change in Asymmetry found by Bendo 
et al. seems to be somewhat in contradiction to the smooth 
decrease found by Lotz et al. The general picture is that as 
long as features smaller than one kiloparsec are resolved in 
the images, the morphological parameters do not suffer too 
much but decline rapidly for coarser images. 

The THINGS observations are designed to sample the 
H I at spatial scales below a kiloparsec but the future all-sky 
surveys will most likely sample the H I disk at lower spatial 
resolutions. To quantify the effects of distance on the mor- 
phological measurements, we use the Hi RO maps of five 
face-on spiral galaxies: NGC 628, NGC 3184, NGC 3351, 
NGC 5457 and NGC 6946. To shift these to a series of dis- 
tances, we cosmologically dimmed the maps, smoothed with 
6" resolution at the appropriate distance, and resampled to 
1" pixel sca le of the THINGS s urvey, adding in the noise 
reported by IWalter et al.1 (|2008D for" a single channel map 
(their table 20. 

In Figure [2] we plot the ratio of the measured morpho- 



^ Our Hi maps are observed throughout the datacube but we 
treat these m aps as single images, sim ilar to the redshifting pre- 
scription from lGiavalisco et al.l 
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logical parameter over the original, fiducial values (RO maps 
at the fiducial distance of zero). 

Concentration and Smoothness of the H I maps are af- 
fected the most by distance effects. In the case of Smooth- 
ness, this is unsurprising as it relies on a smoothed version of 
the original image for the measurement. In the optical, the 
size of the smooting kernel is adapted to be sensitive to spe- 
cific scales of star-formation. Hence, additional smoothing 
with distance is very likely to affect this value. 

Asymmetry, Gini, and A/20 change less than 20% in 
value with increasing distance. This is encouraging for their 
use in throughout a local volume survey. Asymmetry cer- 
tainly does not sh ow the kind of behaviour described by 
iBendio et all l|2007h for the 24 /im SINGS maps (dotted line 
in Figure [2} ■ Considering the 24 /im images trace the dusty 
star-formation and our maps the cool atomic gas, this is not 
as surprising as it initially appears. In the 24 /im maps, a 
large fraction of the flux is in a few star-forming regions, 
dominating morphological measurements. In the case of Hi, 
the morphology is determined by a large area with a very 
limited range of column densities. The difference in contrast 
for these observations may account for their different be- 
haviour with increased distance. EUipticity does not change 
in value much either. These galaxies have very low values of 
ellipticity as they are all near face-on. 

In a companion paper in this series l|Holwerda et al.l 
l201ld ). we quantify the effects of increased distance using 
simulated Hi maps (initially at 150 pc. sampling) and typ- 
ical resolutions, pixels-scales and distances for current and 
planned surveys. From these simulated H I maps, it becomes 
clear that with current and planned observatories, morpho- 
logical Hi measurements are only useful for the very local 
volume. 

3.5 Inclination Effects 

In the case of the Concentration deflned by two circular 
apertures, the inclination dependence is an accepted fea- 
ture. However, even in the case of elliptical apertures, there 
is often a remaining c oncen tration-inclination relation. For 
example iBendo et al.l (|2007l ) find that inclination influences 
the measured Concentration and they introduce a correc- 
tion. Other parameters may well have subtle dependencies 
on the inclination of the disk. To explore these dependencies, 
we use the same five face-on galaxies as for the resolution 
simulations and rotate them around the x-axis assuming an 
infinitely thin disk for the Hi emission (resampling solely 
along the y-axis.). While these galaxies have some intrinsic 
inclination (see Table Al in the Appendix electronic version 
of the manuscript) , we treat them as perfectly face-on as a 
starting point. Figure [3] shows all six parameters as a func- 
tion of inclination. A benefit of Hi observations is a good 
new estimate of the inclination from the tilted ring fit to 
the velocity cube. Any correction for inclination will conse- 
quently be easy to perform for low redshift galaxies. 

Asymmetry, Smoothness and Gini are all relatively un- 
affected below an inclination of 80°, above which our as- 
sumption of a perfectly flat disk becomes tenuous. Ellipticity 
naturally changes drastically by inclining the disk. M20 con- 
verges to a value of -0.7. Concentration f ollows a comp l icated 

pattern peaking in the 30 40° range. iBendo et all (|2007l ') 

flnd a much smoother correction of concentration for incli- 




Distance (Mpc) 



Figure 2. The effect of distance on Hi galaxies. All morpho- 
logical parameter are expresses as tlie ratio between the original 
value and the one found at greater distance. All five galaxies were 
dimmed, smoothed with the resolution of THINGS (6") and re- 
gridded to the THINGS pixel s cale (1" ). The co rrection for Asym- 
metry with distance from Bcndo et al. I h mi l is marked with a 
dotted line (for Spitzer/MIPS data with 6" resolution). 



nation, based on their elliptical apertures. Arguably, they 
pre-corrected for inclination by using an elliptical apertures 
whose shape depends on the inclination of the disk. 



4 RESULTS 

Appendix A {online version of the paper) shows the column 
density maps, plots of the parameters over wavelength and 
lists morphology parameter values for each of our 28 galax- 
ies. The column density map of each galaxy shows the two 
areas in the images, corresponding to the extent of the gas 
(blue) and stellar (red) disk, deflned by us as the 0.3 and 
30 X 10^° cm~^ Hi contours respectively (or alternatively 
0.24 and 24 Mq/pc? respectively). We exclude those parts 
of the images where the H I map is below these values, even 
in the case of Hi "holes" in the Hidisk. All six parameters, 
determined within these contours are plotted for each wave- 
length we obtained data for. In order to compare different 
wavelengths, we deflned the areas over which all the mor- 
phological parameters were computed exclusively based on 
the Hi column density map. Appendix A [online version of 
the paper) also lists our brief notes on what we found for 
each galaxy in the literature, especially regarding its level 
of interaction. First we will compare our parameters to pre- 
vious morphological results and secondly to other measures 
of interaction. 
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Figure 3. The changes in the values of the morphological pa- 
rameters as a function of inclination of the disk. NGC 628, NGC 
3184, NGC 3521, NGC 5457 and NGC 6946 were used. We as- 
sume that the H I disks are at zero inclination at the beginning 
and we rotate around the Y-axes. 

4.1 Comparison to Previous Results 

The SINGS sample, and hence the TH INGS sam p le, ha s 
some overlap w ith the sample from iFrei et ahl (|l99d ). 
IConselic3 (l2003l) pr e sents CAS parameters for the Frei sam- 
ple. and iLotz et al.l (|2004l ') presents Gini, M20 and CAS pa- 
rame ters for the Frei sa mple, both in the optical R-band 
filter. [Bendo et al] (|2007l ) presents Concentration, Asymme- 
try, Gini and M20 for two infrared wavelengths for the entire 
SINGS sample. Here, we compare our results for the H I con- 
tour corresponding to the stellar disk. 

^.1.1 Optical 

IConselic3 (|2003l ) present Concentration, Asymmetry and 
Smoothnes s in t he R-band for the Frei sample, and 
iLotz et ai] (|2004h presents separate estimates of the 
Concentration-Asymmetry-Smoothness values as well as 
Gini and M20 for the Frei sample, in both R and B band. The 
measurements from these two papers are the nearby (iso- 
lated) galaxy re ference fo r CAS measurements. I n Fig ure [4] 
we compare the iConselicd (2003, ) and Lotz et al] |200J) val- 
ues respectively to ours for the R-band values for our "stel- 
lar" disk (the 30 x 10^° cm"^ contour). The main difference 
between our implementation and the earlier work is the def- 
inition of the area over which the parameters are computed; 
we define the area over which the parameters are computed 
based on an Hi contour. In contrast, both previous authors 
use an isophote in the optical or near-infrared imagf^l. The 



Figure 4. The R-band values for CAS from lConselic3 l|2003h (left 
panels) and iLotz et al.l ||2004| ) for the Gini and M20 parameters 
(right panels). 



Frei et al. images and the SINGS R-band ancillary data are 
not exactly the same in depth and filter which may result in 
additional differences. 

There are some small differenc es in the parameters^ 
definitions as well. Both iConselic3 (|2003l ) and iLotz et al.l 
(|2004 ) minimis e the Asymmetry b y all owing the galaxy's 
centre to vary. iBendo et all (j2007l ) and iLotz et all ()2004h 
have a similar a pproach to M2 ■ We chose to use the dynam- 
ical centre from I Walter et al] l|2008l ) and treated variations 
on it as an uncertainty. 

The Gini values from iLotz et aL I l|2004l ) span a much 
smaller range that our values or those from iBendo et al.l 
(2007), possibly because of the difference in area definition. 
Given that our areas are greater than those in Lotz et al., 
there is a relatively bigger contribution by a mass of low- 
intensity pixels, increasing the Gini parameter. Similarly, 
Bendo et al. compute the Gini value for a greater area (the 
RC3 ellipse), including many more low- flux pixels. This dif- 
ference in area would explain higher values of Gini com- 
pared to Lotz et al., but not necessarily why there is a big- 
ger spread in our values. In addition, Bendo et al. note that 
AGN activity change the IR values considerably and exclude 
the centre for some of their galaxies for this very reason. A 
central AGN is likely to be excluded automatically by an H I 
contour. Our concentration values for R-band may well be 
higher is some central pixels are excluded by the H I contour 
(e.g., central Hi hole in NGC 2841, Figure ??). The differ- 
ence in the M20 values is mostly due to the differences in 
area or possibly the depth of the data. 



Optical isophote and surface density contour are technically 
identical terms. We chose to use contour to emphasise that the 



choice of area is done based on the H I column density map, not 
the optical images. 
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4-1.2 Infra-red 

iBendo et alj (|2007l ) present values for Concentration, Asym- 
metry, Gini and M20 for the 3.6 and 24 /im images in the 
SINGS sample and by design there our sample is a subset 
of this sample. , Mu fioz-Mateos et al. (2009) compare these 
infrared values to those in the UV. Bo th authors focus on 
their use for Hubble Type classifiers. iBendo et all (|2007l ) 
correct their Asymmetry values for distance and Concentra- 
tion values for inclination. The other significant difference 
between our determination and that of Bendo et al. is that 
their result is from an elliptica l aperture (defined by the d e 
Vaucouleur's D25 in the RC3 Ide Vaucouleurs et al1ll99lh . 
and ours are for an Hi contour. 

Therefore, we recomputed our implementation of 
the morphological parameters (Concentration, Asymme- 
try aBd_Gmi_and__M2o) for the elliptical apertures used 
in I Bendo et al.l l|2007l ). We pl aced the aper t ures c entered 
on the galaxy position from IWalter et al.l (|2008l '). using 
the minor and major axes from the NASA Extragalac- 
tic Database (http: //ned www. ipac . caltech.edu/ ) and the 



Reference Catalog (RC3 Ide Vaucouleurs et al 



199^ 



posi- 



tion a ngle. Figure [S] shows the values from I Bendo et al.l 
and ours. There is scatter in the values but that is 
to be expected as the input apertures, central positions and 
data are not perfectly identical. Concentration and Asym- 
metry are corrected by Bendo et al. for disk inclination and 
distance respectively. We have uncorrected these values here 
for the comparison. 

Our 3.6 ^m Concentration and Gini values are higher 
than those from Bendo et al. In past, this may be because 
Bendo et al. remove the contribution by the AGN in a few 
of their galaxies, as well as small differences in the adopted 
centre of the disk. There is substantial scatter in the M20 
to higher values for the 3.6 /im and our Asymmetry values, 
in both 3.6 and 24 fim are lower for the highest values for 
Asymmetry from Bendo et al. These differences are because 
we smooth the 3.6 /im data to 6" resolution. In addition, 
some of the differences in M20 and Asymmetry can be the 
result of slight differences on which pixels on the edge of 
the aperture are included in the computation. Both of these 
parameters weigh pixels in favor of those on the edge of the 
object 

All considering, the values for these four parameters 
agree with the previous authors within our computed un- 
certainties. Smoothing effect should be kept in mind when 
comparing this paper's values to others but the smoothing 
is essential for the comparison across wavelength. 



4.2 Morphology of Star-Formation and H I 

Because most of the existing work has been done in wave- 
lengths dominated by recent star formation, mostly UV, we 
compare the morphological parameters determined in the 
Hi RO maps to the parameters we determined in the NUV, 
FUV and MIPS 24 /im. The spatial resolutions in these 
wavelengths are all very similar and we use the same two 
areas; the 3 x 10^^ and 3 x 10^^ cm~^ Hi contours. Figures 
[HI [7] and [S] show the six parameters computed in the far-, 
near-ultraviolet and 24 /im compared to the values in the 
RO Hi maps. 

EUipticity translates extremely well from the UV to the 




Figure 6. Our measurements of Concentration, Asymmetry, 
Smoothness, Gini, M20 and EUipticity in both the stellar disk 
(3 X 10^^ cm~-^, solid points) and the gas disk contour (3 X 10^® 
cm"'^, grey points) in far-ultraviolet compared to the values com- 
puted in the H I RO map. 
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Figure 7. Our measurements of Concentration, Asymmetry, 
Smoothness, Gini, M20 and EUipticity in both the stellar disk 
(3 X 10^^ cm"'^, solid points) and the gas disk contour (3 X lO'^^ 
cm~^, gray points) in near-ultraviolet compared to the values 
computed in the H I RO map. 
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Figure 5. Our measurements of Concentration, Asymmetry, Gini and M20 in an elliptical aperture defined by the major and mino r axes 
(D25) from NED and the RC3 position angle for both the 3.6 and 24 /im images, compared to the values from I Bendo et ahl ||2007| ). The 
dashed line is the line of equality. 




Figure 8. Our measurements of Concentration, Asymmetry, 
Smoothness, Gini, M20 and EUipticity in both the stellar disk 
(3 X 10^^ cm~-^, solid points) and the gas disk contour (3 X 10^^ 
cm~'^, grey points) in Spitzer 24 ^m compared to the values com- 
puted in the H I RO map. 



centre of the galaxies. In the case of the Gini parameter, 
one could expect a shift because the range in values in an 
H I column density map is much smaller than the range in an 
UV flux map. Small changes in Asymmetry and M20 can be 
expected as a bright star-forming region does not necessarily 
translate into a high H I column density. Asymmetry for the 
H I and UV follows each other much better than the 24 /im 
emission. A similar relation is seen in M20 where Hi and 
UV are reasonably similar but the 24 /im is not. Smoothness 
shows no clear relation. 

In general, the morphological parameters from the far- 
or near-ultraviolet translate reasonably well to the H I, with 
the exceptions (Concentration and Gini) well understood. 
The discrepancy is much greater with the 24 pm MIPS 
observations. Thus, Hi and restframe ultraviolet merger 
measurements from morphology translate directly but any 
FIR morphology requirers additional information. Arguably, 
FUV and Hi trace the atomic ISM component and the 24 
fim is closer associated with the molecular {H2) component. 
In favor of this argument is the poor relation in A'ho; the 
bright pixels in 24 ^m do not coincide with high column 
densities in Hi. Fortunately, the majority of morphologi- 
cal measurements of mergers are in rest-frame UV and our 
measurements of merger fractions and rates should translate 
relatively directly. 



5 SIGNATURES OF TIDAL INTERACTION 



Hi. Concentration, M20, and Asymmetry also translate well 
but to a lesser extent. However, this can be expected for 
the Hi disk. For instance, it was already known that gas 
disks are less concentrated, often with depressions in the 



In this section we explore the relation between other meth- 
ods to determine the level of interaction for each galaxy 
to our morphological parameters. The level of interaction 
is difficult to quantify although some authors present a 
parameterisation of the level of tidal force on a galaxy 
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Figure 9. All six morphological parameters, Concentration, 
Asymmetry, Smoothness, Gini, M20 and EUipticity as a f unction 
of tidal disturbance (0) from iKarachentsev et al] ||2004| ) . Grey 
points are for the outer contour, (the gas disk) and the solid 
points for the inner contour (the stellar disk). 



Figure 10. All six morphological parameters. Concentration, 
Asymmetry, Smoothness, Gini, M20 and EUipticity as a function 
of our ranking. Grey points are computed over the gas disk, solid 
points for the inner stellar contour. We note the marked increase 
in Asymmetry. 



(IKarachentsev et al ] |2004l : iBournaud et al]l2005l ). Alterna- 
tively, one can take the level of non-circular motion from the 
THINGS results as a measure of recent interaction. Many 
of these galaxies have been studied extensively and on sev- 
eral, previous authors have remarked signs of disturbance 
(see Table [T]). In this section, we compare our morphological 
parameters in H I to these various parameterisations to iden- 
tify those parameters that appear to be the most promising 
interaction tracers. 

5.1 Tidal Interaction Parameter (O) 

One way to measure the interaction rate is to find close 
galaxy pairs whic h, most lik e ly, wi ll be gravitationally inter- 
acting. Karachentsev et al.l (|2004l ) presents a catalogue of 
neighbouring galaxies and a tidal estimate for each galaxy 
(0). Negative values of O correspond to isolated galaxies, 
and positive values are typical of group members. We should 
note that there are no completely isolated galaxies in our 
sample (9 < 0). Only NGC 628 and NGC 2403 have 6 = 0; 
galaxies in average surroundings and we note that NGC 2403 
was the zero-point calibrator for this tidal parameter. Five 
galaxies have no value for 9 (NGC 2841, NGC 3198, NGC 
3521, NGC 5055, NGC 7331). This makes a direct compari- 
son between the morphological parameters and tidal interac- 
tion more difficult. Figure [9] shows no clear relation between 
any of the morphological parameters of the Hi maps and O, 
compared either over the optical and gas disk. 

Because we only have one galaxy with a high value of 
the Karachentsev tidal estimate (0 = 4 for NGC 5194), 
we can only glean trends with this parameter: compared to 
the locus of group galaxies, the M20 parameter for the H I 



map appears to be higher for NGC 5194, and Concentration, 
Asymmetry and Gini show a slight trend. We interpret the 
relation with O useful to point out those parameters that 
could be of use to identify mergers in Hi morphology, but 
the THINGS sample does not have the spread in O to isolate 
the Hi morphological parameter space where gravitational 
interactions reside. 

5.2 Interaction Ranking 

Because the THINGS sample was mostly chosen to include 
nearby and non-interacting galaxies (M51 being the obvi- 
ous exception), the majority of the above galaxies are not in 
the stage of an interaction where the morph ological signa- 
ture is the strongest (.Lotz et al. ■,2008bl . l2010al lbh. So much so 
that lSmith et al.l (120071 ) use the SINGS sample (minus M51) 
as their reference for non-interacting galaxies. However, we 
can rank the THINGS galaxies on the level of interaction 
signature reported in the literature. In Table [1] we rank the 
sample based on the plausible stage of interaction from iso- 
lated, non-harassed galaxies to merger and merger remnant. 
This ranking is subjective as more information on a galaxy 
(such as high-resolution H I data) often reveals signs of mild 
interaction (e.g., warps or low column density structures), 
and it is still unclear what signatures are from interaction 
or not. For instance, warps may be formed by other means 
as well (e.g., IGM ram pressure). Certain features are long 
lived (e.g., warps) while other fade relatively quickly (tidal 
arms). Figure [10] shows all six parameters as a function of 
the ranking in Table [T] Asymmetry shows a gradient with 
the ranking while M20 does not. In our opinion, this does 
not mean Asymmetry is a better parameter to detect inter- 
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Table 1. The ranking of interaction 







XlOtjGS 




1. 


NGC 2841 


H I warp 




2. 


NGC 3184 






3. 


NGC 3521 






4. 


NGC 3621 






5. 


NGC 3198 


strong non-circular motion 


Trachternach et al. (2008) 


6. 


NGC 2903 


small companion 


Jrwin ct al. f2009) 


7. 


NGC 628 


high-velocity clouds 


Kamphuis & Briggs (1992) 


8. 


NGC 925 


H I tail 


&incisi ct al. (2008) 


9. 


NGC 5457 


lopsided 


Richter & Sancisi ( 1994) 


10. 


NGC 6946 


holes, high-velocity complexes 


Boomsma et al. (2008) 


11. 


NGC4736 


H I streaming in and lopsided. 




12. 


NGC 7793 


Sculptor Group member, warp 




13. 


NGC 7331 


Proximity to Stephans Quintet 


Gutierrez et al. (2002) 


14. 


NGC 2403 


M81 group member. 


Fraternali et al. (2002) 






extra-planar H I 




15. 


Holmberg II 


M81 group member 


Stewart et al. (2000) 


16. 


M81A 


M81 group member 




17. 


DDO 53 


M81 group member 




18. 


Holmberg I 


M81 group member 




19. 


NGC 2976 


M81 group member 




20. 


IC 2574 


M81 group member, 








H I supershell 




21. 


NGC 3031 


M81 central galaxy 




22. 


NGC 3627 


member interacting Leo triplet 




23. 


NGC 5194 


canonical 3:1 interaction 




24. 


NGC 5055 


extended spiral structures. 


Bosma (1981) 






H I warp 




25. 


NGC 4826 


counter-rotating Hi disk 


Braun et aL (1994) 


26. 


NGC 5236 


high- velocity clouds, stellar and 


Malin & Hadlev (1997) 






H I streams and an H I warp. 


Radburn-Smith et al. (2011) 



action, just that it is very sensitive and may well pick up 
other effects or longer lasting effects as well. 



5.3 Non-circular Motions 



iTrachternach et al.l (|2008l ) report a measure of the relative 
non-circular motion in a subset of the THINGS galaxies; 
the total power in the non-circular harmonic components 
over the maximum velocity in the rotation curve {Ar /vmax)- 
They report amplitudes of the non-circular motions lower 
than expected for the steep central mass density predicted 
by ACDM numerical simulations. This THINGS subset of 
19 galaxies was selected to be non-interacting, limiting our 
comparison here. In Figure 1111 we plot the morphological 
parameters as a function of relative non-circular strength. 
At best, there seem to be only very weak trends between 
Concentration, Asymmetry and M20 with the relative non- 
circular motion {Ar). However, the highest value of Ar 
(NGC 3627), is a close group member (there is no Ar value 
for M51). Thus, morphological and dynamical parameters 
combined may delineate a superb space to identify merg- 
ers and harassed galaxies. Future large surveys will have the 
dynamical information automatically assessed by GALAPA- 
GOS opening possibilities in combination with these mor- 
phological parameters. 
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Figure 11. All six morphological parameters. Concentration, 
Asymmetry, Smoothness, Gini, M20 and EUipticity as a function 
of the non-circularity of the H I kinematics in the disk according 
to ITrachternach et al.l 1I2OO8I ). Grey points are the gas disk, solid 
points are computed for the inner stellar contour. 
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6 DISCUSSION 

In this paper we present our determinations of six morpho- 
logical parameters computed within an inner and outer con- 
tour, approximately corresponding to the stellar and gas 
disk of these spirals. We aimed to determine how well H I 
morphology can serve as an indicator of tidal interaction. 
There is much anecdotal evidence for this but the parame- 
terization used to date mostly in restframe ultraviolet im- 
ages offers a way to quantify. We selected the THINGS sam- 
ple as it is nearby, mostly comprised of spirals and offers 
many excellent ancillary data. We tried to answer the fol- 
lowing questions: which disk size (stellar or gas) works best 
to determine the H I morphology in? How well does the H I 
morphology compare to others used for interaction so far? 
Which parameters are the best indicators of interaction in 
Hi? 



6.1 Which wavelength shows interaction? 

From the morphology-wavelength plots in the Appendix A 
{online version of the paper), as well as Figures [6l [T] and 
[HI it is apparent that the wavelength in which the pa- 
rameter is determined is very important. Ultra-violet and 
the Hi parameters appear to correlate as these trace star- 
formation and some of its fuel or alternatively massive stars 
and the atomic co mponent due to their photo-dissociation 
l| Allen et al.ll2004l ). Future Hi results on interaction rates 
of spirals should be very similar to ultra-violet determined 
ones. 

Indeed, the sensitivity of Hi observations to gas-rich 
and minor mergers are a critical benefi t as these types do m- 
inate the mergers at higher redshift l|Lotz et al.l l2010"al lbh. 
The all-sky surveys envisaged, will then provide a valu- 
able local reference for the characterization of these types 
of mergers. The added benefits are that H I appears to be 
very sensitive to gravitational interaction and there is no 
lag time waiting for star-formation to occur. On the other 
hand. Hi morphology may be uniquely sensitive to other 
phenomena, such as ram pressure stripping. 



6.2 Which disk size works best? 

The wavelength in which these parameters are determined is 
critical (see the Figures in Appendix A of individual galaxies 
{online version of the paper) or Figures [HI), but which of 
the two choices to delineate the area for they are computed 
is optimal? There is some extra information in optical mor- 
phological parameters if they are computed over the extent 
of the gas disk (3 x 10^^ cm~'^ contour) as opposed to the 
stellar disk (3 x 10^^ cm~^ contour). Notably, ultra-violet, 
Spitzer 24 fim, and especially of course the H I change with 
increased area. In the case of UV, this is partly because the 
sky level and noise are so low that a signal is easier to pick 
out. The compatibility is good news if we want to compare 
future Hi results, computed of the whole gas disk, to those 
found by other authors on restframe ultraviolet disks. 



6.3 Which parameters are indicators of 
interaction? 

The main limitation to answer this question is the choice of 
sample as the THINGS sample was originally chosen as a 
relatively quiescent reference sample; most THINGS galax- 
ies are group members. With the exception of M51, there are 
no ongoing mergers. In previous studies, both the Asymme- 
try and A/20 parameters were used to identify strongly per- 
turbed systems in HST deep fields. The Hi value of these 
do seem to follow those in the FUV and NUV relatively 
closely, which is promising but not conclusive. If one com- 
pares the tidal parameter or our ranking for the bulk of 
the THINGS sample, there is a gentle rise in Asymmetry 
that one would not expect if this parameter is exclusively 
dependent on interaction yet the A/20 parameter for M51's 
stellar disk does not seem out of place with the rest of the 
THINGS and VIVA[EI sample. Interestingly, the M20 value 
for the entire gas disk does seem to deviate somewhat (Fig- 
ure [§]). Therefore, we speculate here that a combination of 
these two parameters may yield the best result in identify- 
ing merging or strong l y intera cting systems. In our next pa - 
pers, iHolwerda et al.i (|2011bl 'l and iHolwerda et all (|201ld '). 
this speculation is borne out by a large sample of Hi maps 
and simulations of the Hi column density maps during a 
merger. 



7 CONCLUSIONS 

From our multi-wavelength analysis of the quantified mor- 
phology of 28 galaxies in the THINGS survey, we conclude 
the following: 

1. The Hi morphology in column density maps down to 
approximately 10^° cm~^ is promising perspective on galaxy 
interactions. Many of the H I parameters show a close link to 
those in the near- and far-ultraviolet, which have been used 
to date to morphologically identify ongoing mergers. Con- 
centration, Asymmetry and EUipticity appear to be closely 
linked; Smoothness and M20 show a noisier relation between 
H I and UV morphological parameters and the Gini param- 
eter is the notable exception with little correlation between 
the H I and UV determined values. 

2. Despite the close link between far- and near-ultraviolet 
and Hi morphology, the wavelength in which the morpho- 
logical parameter is measured is key, with the choice of area 
over which the parameters are computed of secondary im- 
portance. This has some implications for comparing mor- 
phological parameters across different studies, especially if 
these were done in different wavelengths. The only excep- 
tion seems to be the UV and H I morphologies, which can 
be compared more directly. An increase in sensitivity of the 
H I maps, or a corresponding decrease H I contour level, only 
affects the Gini parameter due to the addition of many more 
low- value data points. The remaining morphological param- 
eters are not affected. This is encouraging for studies com- 
paring Hi morphologies in maps with different depths. 



J^_VLA Imaging of Virgo spirals in Atomic gas, l lChung et al.l 
l200g|) . 
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3. There are, at best, only weak trends between Concen- 
tration, Asymmetry, and M20 and the non-circular motion 
strength {Ar/vmax) (Figure ITlTl . This may imply that Hi 
morphological and kinematic deviations can be complemen- 
tary tracers of disturbances of the spiral disk. 

4. The two most common morphological parameters to 
identify mergers, Asymmetry and M20, show contradic- 
tionary behaviour in Hi. Asymmetry seems to be very sen- 
sitive to disturbances, such as interaction, but also other 
phenomena (starburst etc.), while M20 may be somewhat 
too insensitive (Figures [9l and [TO)) . Based on the THINGS 
sample, which admittedly does not span a sufficient range 
of gravitational interaction, we speculate that a combination 
of these two parameters may be useful for merger selection 
based on Hi morphology. 

In the future, the high-resolution Hi maps of THINGS 
can serve as an excellent reference because these spirals are 
in just the type of small groups that the majority of spirals 
reside in. Bigger Hi surveys such as Westerbork Hi Survey 
Proje ct (WHISP Ivan der Hulst et al.l I2OOII : Ivan der Hulstl 
I2OO2I ). and those to be undertaken with APERTIF, ASKAP, 
EVLA, MeerKAT, and eventually SKA will provide ever in- 
creasing statistics on the Hi morphology of spirals. Using 
some of these bigger samples, we hope to identify the part of 
Hi morphology space which correspond to ongoing mergers. 
In the next papers in this series, we aim to quantify mor- 
phology to the WHISP sample on H I column density maps, 
identify the merger space for H I morphology, the timescale 
mergers are visible, and ultimately infer a merger rate based 
on the WHISP sample. 
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